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BEAM  PROPAGATION  STUDIES  AT  NRL 
JULY  1983  TO  JUNE  1984 
VOLUME  2 


Introduction 

Thin  report  in  comprised  of  simnaries  of  unclassified  papers  presented  at 

the  DARPA/Services  Annual  Propagation  Review,  June  1984,  at  the  Naval 

Postgraduate  School,  Monterey,  California.  A  variety  of  theoretical  and 

experimental  etudlee  are  covered.  Including  a  review  of  propagation  community 

codes  and  topics  in  channel  chemistry,  holeboring,  hose  Instability, 

diagnostics,  and  beam  channel  tracking.  The  papers  and  authors  are  listed  In 

the  order  In  which  they  appeared  at  the  neeting.  This  report  does  not  Include 

those  papers  whose  first  authors  were  Naval  Surface  Weapons  Center  personnel. 
Mssmrrtpf  ippwwwd  July  27.  1984. 


OVERVIEW  OF  BEAM  PROPAGATION  COMPUTER  COOES 
Richard  F.  Hubbard 

JAYCOR,  205  S.  Whiting  St.,  Alexandria,  VA  22304 

Introduction 

Computer  modeling  has  played  a  major  role  in  the  OARPA 
and  Services  beam  propagation  programs.  Beam  stability  and 
transport  involves  coupling  between  complicated  plasma  phys¬ 
ics  and  air  chemistry  processes  which  are  not  readily  treat¬ 
ed  with  purely  analytical  models.  Computer  models  have  been 
used  extensively  to  provide  guidance  in  choosing  design 
parameters  for  present  and  future  experiments.  They  have 
been  particularly  useful  in  suggesting  methods  for  improving 
beam  stability.  This  paper  attempts  to  provide  an  overview 
of  computer  codes  used  in  the  beam  propagation  program. 

Propagation  codes  rely  on  approximate  treatments  of  the 
beam  dynamics,  fields,  and  conductivity  generation.  In  most 
cases,  the  beam  is  assumed  to  be  ultrarelat i vistic  and  par¬ 
axial,  so  vx  <<  v  «c.  The  frozen  approximation,  in  which  z 
and  t  are  replacld  by  z  and  c  -  ct-z,  the  distance  from  the 
beam  head,  is  often  made. 

A  complete  nonlinear  description  in  the  transformed 
variables  would  treat  quantities  as  functions  of  r,  5,  z, 
and  e.  Most  codes  either  assume  3/9  9  -  0  (axisymmetry )  or 
Fourier  analyze  in  9.  In  linearized  codes,  only  one  Fourier 
mode  (usually  m*1)  is  retained;  this  is  appropriate  if  devi¬ 
ations  from  axi-symmetry  are  small.  Channel  tracking  codes 
and  large  chemistry  codes  are  usually  static,  so  3/3 z  -  0. 
Certain  codes  do  not  calculate  radial  variations  self-conslstently 
(1-0  propagation  codes)  or  assiae  3/3;  ■  0  (single  disc  codes). 
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Basic  Components  of  Propagation  Codes 

Most  beam  propagation  codes  contain  modules  for  treating 
beam  dynamics,  conductivity  generation  (air  chemistry),  and 
electromagnetic  fields.  Over  longer  time  scales,  channel 
hydrodynamics  and  radiation  may  also  be  Included. 

Beam  Dynamics:  Axi-symmetric  beam  dynamics  are  usually 
treated  with  particle  simulation,  envelope  equation,  or 
fluid  models.  Particle  simulation  offers  the  most  complete 
treatment  and  does  not  have  to  rely  on  ad  hoc  methods  to 
treat  phase-mix  damping  of  oscillation.  Codes  which  treat 
the  resistive  hose  Instability  add  a  dipole  (m-1)  beam  dy¬ 
namics  model.  The  spread  mass  and  multi-component  models 
partition  eaoh  beam  slice  into  a  number  of  rigidly  oscillat¬ 
ing  components,  while  linearized  particle  simulation  models 
follow  the  trajectories  of  a  large  number  of  particles. 


Field  Equations:  For  models  which  are  not  radially 
resolved,  the  field  equations  are  usually  replaced  by  a 
simple  circuit  model.  The  most  common  radially-resolved 
field  equation  method  uses  the  ultrarelativistic  Lee  field 
equations  and  the  c  and  z  variable  change.  Axl-symmetric , 
linearized  dipole,  and  fully  nonlinear  Lee  field  equation 
solvers  exist  at  a  number  of  laboratories .  The  MRC  field 
solvers  use  the  full  Maxwell  equations  and  do  not  make  the 
ultrarelativistic  approximation. 

Air  Chemistry  and  Conductivity  Generation?  The  simplest 
conductivity  models  assume  constant  momentum  transfer  fre¬ 
quency  v  .  Most  propagation  codes  now  use  more  elaborate 
models  wnich  calculate  a(n  ,  T  )  based  on  rate  equations  for 
n  and  T  or  an  empirica?  T*(E/p)  relationship.  Full  air 
cnemistry  codes  treat  z  as  fixed  and  treat  a  large  number  of 
atomic  and  molecular  species  and  rate  equations.  Nonlocal 
conductivity  codes  allow  beam  generated  plasma  to  be  trans¬ 
ported  and  usually  focus  on  the  role  of  energetic  secondary 
electrons  (delta  rays). 

Channel  Hydrodynamics:  Codes  which  treat  long  time 
scales  must  allow  for  channel  hydrodynamic  motion.  Simple 
adiabatic  expansion  based  on  total  beam  energy  deposition  is 
the  simplest  approach.  A  radially-resolved  1 -D  standard 
hydro  package  is  used  in  some  codes,  while  the  most  elabo¬ 
rate  hydro  treatments  are  the  2-D  and  3“D  turbulence  models 
developed  at  NRL . 

Categorization  of  Computer  Codes 

In  this  section,  various  computer  codes  will  be  grouped 
together  in  broad  categories,  and  a  tabular  comparison  will 
be  made.  Some  codes  may  fall  into  several  categories. 


I .  Axi 

-Symmetric 

Beam 

Propagati 

on  Codes 

Code 

Lab 

Dynamics 

Cond 

Comments 

SIMMO 

Particle 

Fixed  v 

- 

RINGBEARERI 

LLNL 

Particle 

Fixed  v“ 
o(n  ,E/p) 

- 

CPROP 

MRC 

Particle 

Can  treat  low  Y 

SIMPAR 

LJI 

Particle 

o(n?T) 

- 

FLUSH 

SAI 

Fluid 

BMCOND 

like  HIGAP 

FLASH 

SAX 

Fluid 

E/p+LTE 

hydro,  radiation 

FDFC 

LLNL 

Fluid 

Single  disc 

APC-2 

SNL 

Envelope 

BMCOND 

- 

HIGAP 

SAI 

Envelope 

BMCOND 

- 

EPIC/EBFC 

LLNL 

Envelope 

o(n  ,T  > 
BMCOND® 

hydro,  radiation 

SPLASH 

SAI 

Envelope 

shower  physics 

These  codes  can  be  used  to  study  classical  beam  trans¬ 
port,  nose  erosion,  and  in  some  cases,  axi-symmetric  insta¬ 
bilities.  Note  that  moat  dipole  simulations  are  built  upon 
an  axl-symmetr lc  oode  and  can  treat  the  same  problems. 
Fluid  models  of  beam  dynamics  has  received  considerable 
theoretical  intention  in  the  past  year,  but  such  models  have 
not  yet  been  widely  incorporated  into  propagation  codes. 

II.  Linearized  Hose  Instability  Propagation  codes 


Code 

Lab 

Dynamics 

Cond . 

Comments 

VIPER 

NRL 

multi-comp* 

E/p+LTE 

CPT ,  channel  op 

APC-3 

SNL 

multi-comp* 

Fixed  v_ 

new  code 

EMPULSE 

LLNL 

spread  mass* 

«<n  •  t“) 
E/p?LTE 

CPT  option 

PHLAP 

SAI 

spread  mass* 

intense  beams 

SERPENT 

AFWL 

spread  mass* 

o(n  ,  T  ) 

inactive 

SIMM1 

NRL 

particle 

E/p®  ® 

"coarse  grain" 

SIMPAR 

LJI 

particle 

a(n,T) 

orbit  averaging 

RING- 

LLNL 

partiole 

Fixed  \>m 

inactive 

BEARERII 

m 

KM  RAD 

MRC 

particle 

C-indep 

arbitrary  mode 

VALIUM 

LLNL 

Vlasov 

Fixed 

z-indep.  equil. 

•Assumes 

self- 

similar  m*>0  dynamics  (envelope  equation) 

These  have  been  the  "workhorse"  codes  for  the  beam  pro¬ 
pagation  program  since  hose  instability  considerations  often 
dominate  oholoe  of  beam  and  gas  parameters.  Recent  bench¬ 
mark  simulations  among  VIPER,  EMPULSE,  SIMM1  have  given 
excellent  quantitative  agreement  for  lead  pulse  ATA  hose 
predictions  while  providing  insight  into  the  role  of  some 
subtle  conductivity  effects.  High  current  beams  have  also 
received  considerable  attention  with  these  codes  due  to 
their  predicted  favorable  hose  instability  properties. 
Channel  tracking,  magnetically-coupled  multiple  pulse  modes 
and  extremely  long  single  pulses  have  also  been  recently 
studied  with  these  versatile  codes. 

III.  Nonlinear  codes 


mm 


V  *-»  V  TJ. 


0k 


Code 

Lab 

Solver 

Comments 

DYNASTYII 

NRL 

Dynasty 

static;  channel  tracking 

DYNADISC 

SNL 

"Dynasty" 

particle  sim. ;  new  code 

IPROP 

MRC 

CPROP  mod. 

particle  sim.;  few  modes 

LEAD 

NRL 

"Joyce" 

lead  pulse  separation 

- 

LLNL 

"Dynasty" 

faster  field  solver 

’X:Xy>iXy^Ay!y.: 


VS 


The  field  solver  appears  to  be  the  most  difficult  compo¬ 
nent  of  a  nonlinear  code  since  convolutions  involving  many 
Fourier  modes  must  be  evaluated  if  the  deviations  from  axi- 
symmetry  are  large.  The  DYNASTY  field  solver,  which  treats 
the  nonlinear  Lee  equations,  is  probably  too  slow  for  a 
dynamic  code,  so  considerable  effort  has  been  expended  in 
developing  faster  algorithms.  An  algorithm  developed  by 
Joyce  looks  particularly  promising  but  is  still  being 
tested . 


IV.  Channel  Hydrodynamics  Codes 


Code 

Lab 

Dimensionality 

Comments 

EPIC/EBFC 

LLNL 

r , 

C,  z 

Envelope  beam 

PHAZR 

NRL 

C. 

limited  z 

Model  turbulence 

FLUSH 

SAI 

r » 

C»  z 

fluid  beam 

HINT 

NRL 

r. 

C 

CHMAIR  interface 

FAST  2D(3D) 

NRL 

*. 

y,  c,  (z) 

Turbulence  studies 

HYDRO 

LJI 

r  • 

c 

SIMPAR  interface 

Hydrodynamics  codes  treat  channel  motion  on  long  time 
scales.  EPIC  and  FLUSH  have  been  primarily  used  to  treat 
long  CPT  pulse  trains  and  thus  must  include  radiation  trans¬ 
port  and  high  temperature  chemistry.  The  other  codes  have 
been  designed  to  study  the  formation  of  NIPS  channels.  Note 
that  c  is  the  time-like  variable  in  these  codes. 

V.  Beam  Conditioning/IFR  Codes 


Code 

Lab 

Beam  Dynamics 

Comments 

BEAMFIRE 

LLNL 

Particle 

(m-0) 

Treats 

IFR,  foil,  and 

air 

CENTRING 

LLNL 

Particle 

(m-1 ) 

Simpler 

chemistry 

"Hire" 

LLNL 

Particle 

Single 

disc;  external 

B 

"IFR" 

SAI 

Envelope 

♦ 

For  IFR 

track ing 

rigid  displace 

These  codes  are  designed  primarily  for  wire  zone  and  IFR 
conditioning  cell  design.  BEAMFIRE  predicts  the  beam  radius 
and  emlttance  of  a  beam  as  it  emerges  into  the  air  from  an 
entrance  foil,  while  the  other  codes  primarily  treat  elec¬ 
trostatic  tracking  and  damping  of  transverse  oscillations. 


VI.  Nonlocal  Conductivity/Delta  Ray  Codes 


I 


Code 

Lab 

Method 

Comments 

v.v.v 

•V-Vv 

Nuts 

LLNL 

Boltzmann 

Relativistic  Upgrade 

V-.‘ 

*  •/  %  * 

CAS1 

SAI 

Monte  Carlo 

Emphasizes  6*s 

LOCHEM 

Ski 

Multi-fluid; 

Ageing  model  j's 

• . 

CPROP 

MRC 

particle  sim. 

Also  follows  beam 

*  •  * 

Nonlocal 

conductivity  processes 

can  play  a  major  role, 

especially  at 

reduced  pressures. 

Current  multiplication  in 

Code 

Lab 

CHMAIl;  II 

NRL 

HICHEM 

SAI 

HITECH 

NRL 

ETA  appears  to  Involve  these  effects. 


VII.  Full  Air  Chemistry  Codes 

Comments 

Some  versions  radially-resolved; 
adiabatic  hydro  option 
detail  balance,  high  T 

High  T  version  of  CHMXIR 
S 

This  is  an  abbreviated  list  of  the  full  scale  air  chem¬ 
istry  codes  discussed  in  the  review  by  All.  Such  codes  are 
similar  in  scope  to  the  nonlocal  chemistry  codes  but  provide 
much  more  detailed  treatment  of  the  various  air  chemistry 
reactions  and  species. 

Summary 

Several  important  trends  have  emerged  in  the  past  year 
which  should  indicate  the  planned  emphasis  of  code  develop¬ 
ment  in  the  next  year.  Nonlinear  codes  are  currently  under 
development  at  LLNL ,  NRL ,  SNL ,  and  MRC  and  are  likely  to 
absorb  much  attention  in  the  future.  Magnetically-coupled 
propagation  modes  such  as  CPT  are  just  beginning  to  receive 
detailed  attention,  and  the  complicated  chemistry  and  radia¬ 
tion  processes  which  accompany  such  intense  beams  will  re¬ 
quire  considerable  further  effort.  Fluid  beam  codes  show 
promise  and  may  also  play  a  major  role  in  the  future.  Few 
existing  codes  can  treat  the  hollow,  rotating  beam  which 
emerges  from  RADLAC,  but  the  theoretical  foundations  for 
such  codes  have  been  developed  considerably  during  the  past 
year.  Finally,  the  approaching  propagation  experiments  on 
ATA  and  RADLAC  should  provide  the  first  real  test  of  the 
ultrarelativistic  codes,  and  such  codes  should  play  a  major 
role  in  interpreting  experimental  data  and  suggesting  exper¬ 
imental  parameters. 
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BEAM  IHDUCED  COHDUCTIVITr  EXPERIMEHTS 

M.  Raleigh,  D.P.  Morphy,  R.E.  Pechacek,  and  J.R.  Greig 
Baral  Research  laboratory 
Plasm  Physics  Dirts  ion 
Experimental  Plasm  Physics  Branch 
Washington,  DC  20375 

arraoDocnoH 

An  slsctron  beam  ionizes  the  air  through  which  it  passes.  This  beam 
induced  conductivity  in  turn  det ermines  how  the  electron  beam  propagates.  How 
subsequent  pulses  will  propagate  is  determined  by  the  decay  of  this 
conductivity. 

In  this  experiment  we  hare  looked  for  the  effect  of  water  vapor  on 
conductivity  at  early  times  (<  100  ns)  by  observing  the  influence  of  water 
content  on  the  net  current  (1^  +  I.  A  gnnt ) .  These  measurements  have  been 

performed  in  air  and  nitrogen  with  eathgaa  being  tested  both  dry  and  wet, 
l.e.  •  with  one  molecule  of  water  to  each  55  +  5  molecules  of  air  or 
nitrogen.  Pressures  of  one  atmosphere  and  kO  Tbrr  have  been  tested.  At  one 
atmosphere,  one  part  In  55  represents  -  60%  relative  humidity.  This  ratio  at 
hO  Tbrr  represents  the  conditions  that  might  exist  subsequent  to  the 
hydrodynamic  expansion  of  a  beam  heated  channel. 

APPARATU3 

The  beam  source  is  a  Pulserad  310  generator  driving  a  cold  cathode 
diode.  The  diode  is  coupled  to  the  drift  tube  via  an  «„m  transport  tube 
(length  50  cm).  Using  the  transport  tube  a  reproducible  13.5  kA,  1  MeV,  22  ns 
FVOf  beam  was  injected  into  the  drift  tube  through  a  2.3  cm  aperture  (Figure 
1).  The  transport  tube  was  always  filled  with  20  Tbrr  of  dry  !«.  A  small 
Pbgovskl  coll  was  mounted  Just  upstream  of  the  titanium  foil  separating  the 
transport  tube  from  the  drift  tube.  This  coll  monitored  the  shot  to  shot 
reproducibility  of  the  injected  net  current. 

The  drift  tube  in  which  the  net  current  measurements  were  made  was  an 
aluminum  cylinder  60  cm  long  by  60  cm  in  diameter.  In  a  preliminary  set  of 
experiments  it  was  determined  that  the  net  current  would  be  independent  of  the 
drift  tube  length  for  lengths  over  ~  25  cm  and  that  the  net  current  decreased 
slightly  with  Increasing  drift  tube  radius  for  radii  over  ~  25  cm.  (This  is 
to  be  expected  as  a  result  of  the  increasing  return  path  inductance.)  The 
drift  tube  used  was  choosen  large  enough  to  avoid  the  anomalies  associated 
with  the  small  dimensions  but  small  enough  that  a  lumped  Inductance  would  be  a 
good  model  for  the  drift  tube  as  a  current  return  path.  (Double  transit  time 
in  the  tube  •  4  ns,  net  current  risetime  5  12  ns.) 

Fof  the  net  current  measurements  presented  here  a  set  of  eight  equally 
spaced  B  loops  were  Installed  in  a  shielded  groove  in  the  end  plate  of  the 
drift  tube.  These  loops  were  summed  in  parallel  to  yield  a  fast 
(t  4  «  300  ps)  Rcgowski  belt  of  25  cm  diameter.  (A  60  cm  diameter  Rogovski 

belt  nas  subsequently  been  installed  and  is  described  elsewhere  in  the 


proceedings  of  this  meeting.  The  25  cm  diameter  belt  operates  on  similar 
principles.)  Tb  both  check  the  calibration  of  the  Rogovski  coil  and  measure 
the  injected  HEB  current  (Figure  2),  measurements  were  made  using  a  very  short 
drift  tube  (~  3  cm)  terminated  with  a  Faraday  cup  and  evacuated  to  <  1  a  Tbrr. 

To  supply  the  various  water  /gas  mixtures  to  the  drift  tube  a  flow-through 
system  was  adopted.  The  gas  (water  punped  Hg  or  zero  grade  air)  was  supplied 
from  a  regulated  cylinder.  When  wet  gas  was  desired  it  was  passed  through  a 
water  saturated  sponge  inside  a  15  cm  diameter  lucite  tube.  Ihe  flow  rate  and 
pressure  in  the  lucite  tube  remained  constant  regardless  of  the  particular 
experiment  being  performed.  The  gas  supply  was  throttled  to  the  desired 
pressure  as  it  entered  the  drift  tube.  Gas  was  continuously  pumped  from  the 
tube  through  another  throttle  valve.  Bulk  flow  was  maintained  through  the 
valves  so  that  the  cosposition  of  the  mixture  was  not  altered.  The  water 
content  was  determined  by  observing  the  temperature  at  which  water  or  ice 
appeared  on  a  chilled  mirror  surface. 

RESULTS 

All  net  current  measurements  in  the  drift  tube  have  been  normalized 
against  the  average  injected  net  current.  Thus  the  traces  shown  in  Figures  3 
and  4  represent  average  net  currents  for  each  condition  and  may  be  conpared 
directly  with  one  another.  Four  or  five  shots  were  averaged  numerically  for 
each  condition  and  the  doted  lines  show  the  RMS  spread  in  the  data.  Since  the 
difference  between  wet  and  dry  was  the  same  for  nitrogen  as  for  air  only  the 
'results  for  air  are  shown. 

The  following  conclusions  may  be  drawn  from  a  comparison  of  the  average 
net  current  traces: 

1.  At  40  Tbrr  in  either  ^  or  air  there  is  no  difference  between  the  gas 
being  wet  or  dry  for  the  net  current  during  the  rising  portion  of  the 

pulse  (0  ♦  30  ns).  There  is  a  slightly  higher  net  current  in  wet  gas  during 
the  tail  of  the'pulse ,  but  even  this  difference  is  less  than  the  RMS  spread  in 
the  measurements. 

2.  At  one  atmosphere  in  either  N2  or  air  net  current  peak  is  higher 
and  occurs  later  in  dry  gas.  After  the  current  peak  the  net  current  is  always 
higher  in  the  dry  gas. 

3*  The  net  current  is  slightly  less  in  air  than  in  N2  under  all 
conditions. 

4.  For  either  gas  the  peak  net  current  is  significantly  less  in  40  Tbrr 
than  in  1  atmosphere  and  the  decay  time  is  ouch  longer. 
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On  the  Vibrational  Relaxation  and  Dissociation  of  Nitrogen 
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1.  Introduction.  There  are  several  reasons  why  an 
understanding  of  the  vibrational  excitation  processes  in 
nitrogen  is  lnportant  for  the  propagation  of  an  intense 
electron  bean  through  the  air.  On  the  one  hand,  the  large 
oross  sections  for  vibrational  excitations  with  plasma 
electrons  aeans  that  an  appreciable  amount  of  the  energy 
deposited  by  the  beam  ends  up  in  vibrations.  This  keeps  the 
eleotrons  cooler  which  inhibits  ionization  and  enhances 
recombination.  On  the  other  hand,  the  relaxation  of  this 
energy  into  the  kinetic  mode  is  relatively  slow.  The  rate 
of  release  is  critical  in  determining  channel  formation 
characteristics.  Furthermore,  the  role  of  electrons  in 
dissociating  the  molecules  from  higher  vibrational  levels 
and  the  complete  dissociation  of  N.  in  high  current  modes 
requires  a  detailed  study  of  the  vibrational  energy  storage 
and  its  relaxation. 

2.  Models  and  Results.  To  illustrate  some  of  these  ideas 
we  will  give  some  typical  results  from  runs  of  the  NRL 
CHMAIR-II  code  and  briefly  describe  its  vibrational  energy 
model.  Then  we  will  dlsouss  a  detailed  nitrogen  chemistry 
code  VBRAD-II  which  was  developed  specifically  to  address 
the  problems  mentioned  above. 

Tables  la  and  1b  show  the  results  in  one  atmosphere  of  N_ 
uslng  CHMAIR-II.  A  beam  with  current  of  10  kA,  a  rise  ana 
fall  time  of  1  ns  each  and  a  plateau  of  8  ns  was  utilized 
for  two  oases  of  beam  radius  R^  •  2  sm  and  5  mm. 

Table  la  shows  several  quantities  at  5  ns,  half  way  through 
the  pulse.  The  eleotron  temperature  is  .higher  than  the 
vibrational  temperature,  while  the  gas  temperature  is  much 
lower.  Roughly  tog  of  the  deposited  energy  is  in 
vibrational  exoitatlon.  Table  1b  shows  the  same  quantities 
at  30  ns.  By  this  ‘  time  the  electron  and  vibrational 
temperatures  are  nearly  equal  and  around  60£  of  the  energy 
is  held  in  the  vibrational  mode. 
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Table  la.  5  na 
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D 


V  v. 


V5 


N 

CHMAIR-II 

VERAD 

CHMAIR-II 
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VIB.TEMP(eV) 

.31 

.29 

.11 

.11 

GAS  TEMP (eV) 

.042 

.055 
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.028 
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ELEC  TEMP  (eV) 

.45 

.45 

.32 
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N-DENSITY 

4.8(17) 

5.0(17) 

6.8(16) 

7.0(16) 
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4.5(16) 
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1.1(16) 

1.1(16) 
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49 

37 
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13 
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Table  1b.  30  aa 
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VIB  TEMP(eV) 
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.65 

.18 

.16 
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GAS  TEMP(eV) 

.065 

.090 
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.034 
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ELEC  TEMP(eV) 
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.65 

.23 
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65 

61 

55 
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6 

12 

6 

13 
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2.1  CHMAIR-II  Model.  The  vibrational  aodel  used  in  CHMAIR- 
II  la  baaed  on  the  assumption  m_  is  a  harmonic  oscillator 
and  that  the  vibrational  levels  are  in  a  Boltzmann 
distribution  characterized  by  a  vibrational  temperature  T 
due  to  the  rapid  exaltations  and. de-excitations  by  electron 
impact.  The  exitation  rates i  are  obtained  using  the 
measured  cross  sections.  Thus  the  vibrational  energy  per 


moleoule , 
<t  evM2 
dt  * 


can  be  written 
8 


•oMeN2 


iX1(Exp(le0(T#-T7)/TtfTy)-1) 
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where  X.  »  X,(T  )  la  the  rate  for  electron  impact  excitation 
of  1  quanta,  e  *»0. 3  eV  la  the  vibrational  energy  apaclng  and 
T  and  T  are  the  electron  and  vibrational  temperatures .  N 
la  the  electron  density  and  Ng  is  the  nitrogen  density.  * 

The  vibrational  relaxation  due  to  heavy  particle  collisions, 
e.g.,  N2  0,,  and  0  ap«  added  to  Eq.  1  independent  of 
electron  effects  due  to  different  time  scales  in  their 
collision  rates. 


2.2  VERAD-I  and  -II  Models.  To  ascertain  the  assumption  of 
the  Boltzmann  distribution  for  vibrations  and  to  include  the 
effeots  of  heavy  partiole  collisions  on  the  vibrational 
levels  we  developed  VERAD-I  and  -II  to  describe  the 
vibrational  energy  and  its  relaxation  using  a  master 
equation  that  solves  for  32  levels.  Three  types  of 
Interaction  are  considered  in  this  model:  VE,  VT  and  VV. 
The  VE  interactions  are  the  excitation  and  de-excitation  of 
vibrational  levels  by  electron  impact .  The  VT,  vibrational- 
translational  ,  interactions  are  those  in  which  a  vibrational 
quantum  is  gained  or  lost  during  a  collision  with  a  heavy 
particle.  The  VV,  vibrational-vibrational,  interactions  are 
collisions  between  nitrogen  molecules  in  which  a  quantum  is 
exchanged.  The  details  of  the  master  equation  for  VERAD-1 
and  the  appropriate  rates  can  be  found  in  Reference  5. 

In  VERAO-I  we  assumed  a  constant  electron  density  and 
constant  (but  not  equal)  electron  and  gas  temperatures. 
This  model  shows  the  development  of  the  vibrational  level 
distribution  which  is  close  to  Boltzmann  and,  hence, 
verifies  the  CHMAIR-II  model  for  the  vibrational 
temperature . 

VERAD-1I  solves  the  master  equation  including  the  electron 
beam  interaction  with  the  vibrational  levels  (including 
ionization  and  dissociation).  Furthermore,  the  electron  and 
gas  temperatures  are  calculated  from  first  principles.  The 
ohmic  heating  is  obtained  using  a  circuit  equation. 

3.  Discussion  of  Results.  Tables  la  and  b  give  the  VERAD 
results  for  the  examples  discussed  earlier.  The  vibrational 
temperature  calculated  by  CHMAIRs-II  compares  quite  well 
with  that  of  VERAO.  VERAD  predicts  slightly  more 
dissociation.  This  is  oonsistant  with  the  enhanced 
dissociation  rates  of  vibratlonally  exoited  molecules.  The 
gas  temperature  predicted  by  VERAD  is  higher.  The  reason  is 
that  CHMAIR  II  oarrles  several  of  the  electronically  excited 
states  of  N2  *nd  H  whereas  VERAD  assumes  that  energy  is 


Instantaneously  quenched.  The  electron  density  in  CKMAIR  is 
higher  at  later  times  because  it  uses  a  more  sophisticated 
colllsional<-radiati  ve  recombination  model. 

finally  we  show  one  example  to  illustarte  the  processes 
oontributing  to  dissociation  and  the  effect  of  including 
quenching  due  to  atomic  oxygen.  Figure  la  shows  the 
electron,  vibrational  and  gas  temperatures  of  30  kA,  5  mm 
beam  with  a  pulse  length  of  50  ns  in  pure  nitrogen.  Though 
the  vibrational  and  electron  temperatures  equilibriate  soon 
after  the  pulse  passes,  the  gas  temperature  is  considerably 
lower  at  a  microsecond.  If  we  include  a  fixed  amount  of  3  x 
1  o'  per  oc  of  oxygen  atoms,  about  the  amount  predicted  by 
CHMAXR  for  the  same  beam  in  air,  equilibrium  is  nearly 
reached  at  one  microsecond  as  shown  in  Pig.  lb.  At  the  end 
of  the  run  in  Fig.  la  the  following  processes  contributions 
to  dissociation  were: 

direct  production  of  N  by  the  beam 
production  H  by  dissociative 
recombination 

dissociation  by  secondary  eleotron 
impact ,  and 

dissociation  by  exaltation  to 
32nd  level 
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Introduction 

The  stable  propagation  of  a  aultlpulse  beam  over  long  distances  In  a 
gas  depends  heavily  on  the  hot  channel  which  the  bean  fonts.  If  the  channel 
quickly  reaches  a  proper  balance  between  gas  density  and  conductivity  over 
the  first  few  pulses,  the  reaalnlng  pulses  will  be  able  to  reach  an 
equilibrium  condition  while  experiencing  a  minimum  of  expansion  due  to 
Mordsleck  scattering.  After  a  bean  pulse  has  traveled  through  a  given 
region,  a  significant  fraction  of  the  energy  which  has  been  deposited  in  the 
gas  In  that  region  will  reside  In  molecular  vibrational  states.  To  properly 
model  the  physics  of  beam  channels,  a  detailed  treatment  of  nonequilibrium 
air  chemistry  is  therefore  necessary.  Toward  that  end,  All  et  al.  have 
developed  computer  codes  for  both  wet  (1]  and  dry  air  chemistry  [2]  and  have 
studied  the  interactions  among  the  various  channel  and  beam  constituents 
extensively  [3].  Recently  All  and  Sllnker  have  developed  a  radially 
resolved  (one -dimensional)  version  of  the  CHMAIR  II  code  in  order  to  model 
channel  evolution  within  a  few  hundred  nanoseconds  after  a  given  pulse  has 
passed  a  point  along  the  beam  trajectory. 

Also  of  great  Importance  In  modeling  the  development  of  the  beam 
channel  Is  the  hydrodynamic  expansion  of  the  hot  gas  over  a  timescale  of 
microseconds  following  a  beam  puls*.  Such  expension  cools  the  channel 
between  pulses  and  thus  can  significantly  affect  the  composition  of  the 
channel  gas.  To  account  for  this.  All  and  Sllnker  incorporated  a  fluid 
model  besed  on  adiabatic  expansion  Into  the  zero-dimensional  CHMAIR  II  code 
and  computed  channel  properties  at  the  axis  for  a  WIPS  pulse  train  [4].  In 
this  paper,  we  report  the  development  of  a  radially  resolved  Interpulse 
air  chemistry  model  In  which  fluid-dynamic  effects  are  computed  using  the 
more  rigorous  method  of  flux-corrected  transport.  After  describing  the  code 
in  some  detail,  we  present  a  calculation  of  channel  properties  for  a  pulse 
train  with  typical  ATA  parameters.  This  code  should  be  quite  useful  in 
analyzing  ATA  and  RADLAC  experiments  and  the  WIPS  parameter  space. 

Interpulse  Chemistry  Model 

Our  new  model  computes  the  radial  distribution  of  channel  properties  in 
a  given  plane  transverse  to  the  trajectory  of  the  beam.  The  unique  feature 
of  the  code  Is  the  merging  of  a  state-of-the-art  compressible  flow  algorithm 
(flux-corrected  transport  [5])  with  an  appropriate  chemical  reaction  scheme 
In  order  to  compute  the  evolution  and  composition  of  the  channel  in  detail. 
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For  this  reason,  vs  call  ths  cods  "HINT",  which  stsnds  for  "hydrodynamics 
sad  intsrpulso  chemistry. "  Ths  calculation  proceeds  by  the  method  of  time- 
step  splitting,  in  which  the  equations  for  conservation  of  mass,  momentum, 
total  energy,  vibrational  energy,  and  species  number  are  first  Integrated  by 
flux-corrected  transport  for  s  given  time  step  with  the  chemistry  frozen. 

Ths  resulting  values  then  serve  as  initial  conditions  for  an  implicit 
integration  of  the  chemical  rate  equations  over  the  same  times tep  while 
the  fluid  motion  is  frozen.  We  currently  assume  charge  conservation  in  each 
grid  cell  in  order  to  compute  the  electron  density,  and  we  set  the  electron 
and  vibrational  temperatures  equal.  Detailed  calculations  using  CHMAIR  II 
have  verified  the  validity  of  these  assumptions.  Currently  the  code  does 
not  account  for  beam-generated  turbulence. 

To  integrate  the  chemical  rate  equations,  we  use  a  vectorized  version 
of  CHEHEQ  [6],  which  is  especially  fast  snd  accurate  in  the  nonequilibrium 
situations  characteristic  of  beam  channels.  In  addition,  we  use  an  auto¬ 
matic  rate  processor  to  construct  the  chemical  rate  equations  from  literal 
data,  which  depict  the  actual  chemical  reactions,  and  numerical  data  which 
describe  the  rate  coefficients.  Thus,  changing  the  chemical  reaction  scheme 
requires  no  programming,  unlike  the  usual  method  of  hardwiring  a  specific 
set  of  reactions  into  the  chemistry  package.  In  fact,  a  recent  version  of 
the  rate  processor  permits  one  to  enter  the  temperature  equations  into  the 
code  as  a  set  of  pseudoreactions.  The  various  temperatures  then  become 
additional  species  in  the  reaction  scheme. 

The  Interpulse  chemistry  model  currently  includes  the  following 
variables:  gas  temperature,  vibrational  energy  per  nitrogen  molecule,  and 
number  densities  for  M,  0,  02***  and  ®*  The  Interpulse 

chemistry  currently  consists  of  24  chemical  reactions  of  ths  following 
types:  neutral  association,  rearrangement,  dissociative  recombination, 
neutral  ionisation,  ion-ion  neutralisation,  attachment,  and  detachment. 

Initial  Calculation 

To  test  and  calibrate  HIST,  we  have  considered  a  train  of  five  ATA 
pulses  with  the  following  parameters:  (1)  pulse  radius — 5  am,  (2)  energy — SO 
MaV,  (3)  beam  current — 10  kA,  (4)  rise  time — 1  ns,  (5)  pulse  length—10  ns, 
and  (6)  pulse  repetition  lntervsl— 30  us.  The  calculation  begins  with  a 
simulation  of  the  first  100  ns  by  radially  resolved  CHMAIR  II.  This  time 
interval  is  short  with  respect  to  the  hydrodynamic  time  scale  r/c.-lO  us, 
so  that  we  would  expect  mialaal  error  from  the  exclusion  of  hydrodynamic 
motion  by  that  code.  The  channel  conditions  at  100  ns  serve  as  initial 
conditions  for  ths  HIM!  code,  which  simulates  the  remainder  of  the  pulse 
repetition  interval.  The  channel  conditions  at  that  time  become  the  initial 
parameters  for  the  CHMAIR  II  calculstlon  of  energy  deposition  by  the  next 
pulse.  In  this  manner,  we  attempt  to  give  the  most  accurate  treatment  of 
the  two  stages  of  holeboring:  energy  deposition  and  hydrodynamic  expansion. 
The  results  correspond  to  a  position  near  the  nozzle  from  which  the  beam 
exits. 

Figure  1  shows  aass  density  versus  radius  30  us  after  the  fifth  pulse. 
The  shock  waves  from  this  and  the  previous  two  pulses  are  visible  at  radii 


of  approximately  2.0,  3.0,  and  4.0  ca  from  the  center  of  the  channel,  and 
the  minimum  density  la  roughly  10  per  cent  of  the  ambient  value, 
corresponding  to  a  Maximum  channel  temperature  of  0.32  eV.  Ve  remark  that 
the  channel  la  not  yet  at  pressure  equilibrium  with  the  ambient  air  after  30 
vm;  in  fact  an  average  overpressure  of  0.3  atm  still  exists.  Thus  models 
which  assume  adiabatic  expansion  should  use  an  equilibration  time  scale  of 
about  60  m  for  the  later  pulses. 

Figure  2  shows  the  value  of  the  heavy  particle  density  at  the  channel 
axis  30  m  after  each  pulse  and  the  line  gives  corresponding  values  from  a 
recent  calculation  by  SI inker  and  All,  in  which  the  gas  temperature  dropped 
according  to  the  adiabatic  gaa  law.  The  latter  used  CHMAIR  II  to  compute 
the  value  of  the  density  only  on  axis  and  did  not  provide  radially  resolved 
profiles.  The  agreement  between  the  two  is  remarkable,  although  the  HINT 
channel  is  slightly  cooler,  since  some  of  the  energy  goes  into  shocks. 

Figure  3  shows  plots  of  conductivity  on  axis  versus  time  after  energy 
deposition  for  pulses  3,  4,  and  3.  The  aggreenent  with  the  calculations  of 
Slinker  and  All  la  good  for  pulses  3  and  4  while  the  axis  conductivity 
computed  by  HINT  is  approximately  SO  per  cent  higher  for  pulse  5.  We  are 
still  investigating  this  discrepancy  at  tbs  present  time. 

The  results  of  this  calculation  Indicate  that  the  HINT  pode  should  be 
quite  useful  in  analysing  channel  chemistry  end  the  holeboring  properties  of 
particle  beams  propagating  in  gaaea.  We  will,  however,  require  further 
testing  of  the  convergence  properties  of  the  code  under  various  conditions . 
We  have  found,  for  example,  that  the  code  la  quite  sensitive  to  the  manner 
in  which  the  vibrational  energy  is  transported  and  that  noise  in  the 
pressure  la  enhanced  by  the  chemistry  portion  of  a  tlmestep.  These 
properties  appear  to  be  more  important  for  later  pulses,  and  a  thorough 
study  will  be  necessary  to  understand  the  effects  of  these  features  on  our 
results. 
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Effects  of  Temperature-Dependent  Conductivity  on  Hosa 


S.  Slinker*,  R.  Hubbard*,  R.  Fernaler, 
A.  All,  and  M.  Lampe 


Naval  Research  Laboratory 
Washington,  DC 


1 .  Introduction.  In  this  paper  we  examine  the  effects  of 
dipole  air  chemistry  on  the  resistive  hose  instability.  The 
main  emphasis  will  be  on  the  variation  of  the  electron 
neutral  collision  frequency  v  with  respect  to  the  electron 
temperature.  This  seems  to  be  the  most  important  "exotic" 
air  chemistry  effect  in  the  nose  of  ATA-type  beams  in  full 
density  air.  For  these  beams  there  is  little  avalanche 
ionization.  The  hose  instability  is  most  dangerous  in  the 
nose  region  where  the  conductivity  is  low.  It  will  be  shown 
that  dipole  chemistry  effects  act  to  reduce  the  effective 
conductivity  in  this  region,  thereby  allowing  the  hose  to 
grow  more  than  if  these  effects  were  neglected. 


2.  Collision  Frequencies.  If  q^v)  is  the  cross  section  for 
momentum  transfer  for  collisions  of  electrons  of  speed  v 
with  a  species  of  particle  density  M  and  f(v)  is  the 
velocity  distribution  function  then  the  collision  frequency, 
v,  is  given  by 


1  - 


v<jM  by 


where  n  is  the  electron  density.  If  q,  ls^gnstant  and  f  is 
Maxwellian  of  temperature  T  then  v  “  M?  _3/_When  M  is  an 
ion  the  Coulomb  interaction  gives  v  *  MT-'3  if  the  weak 
dependence  of  InA  on  T  and  M  is  ignored.  In  Fig.  1  we  show 
the  collision  frequency  coefficients  for  nitrogen,  oxygen 
and  air  calculated  from  cross  sections  compiled  in  Ref.  1. 
The  typical  electron  temperature  in  a  beam  is  on  the  order 
of  1  eV .  For  moderate  values  of  the  eletric  field  (we 
ignore  magnetic  and  ac  effects)  the  conductivity  is  given  by 

2  _ 


tant 


.-3/2 


When 


is  an 


In  discharge  studies  where  there  is  a  constant  electric 
field,  E,  the  concept  of  a  drift  velocity  v  «v.(E)  is  often 
used.  In  this  ease  the  plasma  current  density  is  given  by 
J  »env.«cE  so  c-env.(£)/E.  Although  the  induced  electric 
field  Is  not  constant  in  Intense  beams,  its  variation,  for 


J  •env.-cE  so  ff-env.(£)/E. 
field  Is  not  constant  in 
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slow  rise  times,  is  on  the  order  of  a  nanosecond  while  the 
collision  frequency  for  the  plasma  electrons  is  typically 
10T  in  full  density  air. 


3.  Effects  on  Hose.  The  collision  frequency  depends  on  the 
eleotron  temperature,  which,  in  turn,  depends  strongly  on 
the  magnitude,  E,  of  the  electric  field  because  of  ohmic 
heating.  In  the  following  we  assume  that  the  collision 
frequency  is  a  function  of  E  explicitly. 

The  beam  body  is  the  region  behind  the  pinch  point  but 
before  the  region  where  recombination  processes  begin  to 
become  important.  In  the  beam  body  Maxwell's  equations 
reduce  to  Ampere's  law* 


-»?*  •  H  <jb  *  y 

Where  A  is  the  axial  component  of  the  vector  potential  J.  ii 
the  beam  current  density  and  J  »aE  •-e^A/dc  is  the  plasma 
current  density.  g«v t-z  is  the  distance  behind  the  beam 
head .  z 

The  perturbation  of  this  equation  leads  to  the  following 
equation  for  the  dipole  quantities. 


b  p  br*  -  *  H  tJ< 


♦  v 


where  dipole  quantities  are  denoted  by  .  Using  v-w(E)  with 
S-lEJ  the  magnitude  of  the  electric  field  in  (1)  the 
perturbed  plasma  current  becomes,  after  some  algebra. 


-  aE. 


ffE  •  c^E  -  ®{ 1 

°n6z  \»  dE1 


Substituting  this  in  (3)  gives  .  A 

^  t  a  dA  .  _  n. 

5?  7  fr  rA  '  ~  s.ff  a?  •  *  ~'Jb  *  ,E,  «>  (5) 

where  the  effective  conductivity  is  defined  by 

*.rf  -  *1'  -  *  It1-  (5) 

In  terms  of  the  drift  velocity  o#ff«en  dv^/dE. 

In  full  density  air  for  the  nose  of  ATA-type  beams  dv/dE  is 

positive  (Fig.  1),  so  the  effective  conductivity  is  lower 

than  what  it  would  be  if  v  were  constant.  This  leads  to  an 


expectation  of  more  hose 
supported  by  code  results 
seotion. 


growth  and  this  conclusion  is 
from  VIPER  given  in  the  final 
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If  v  varies  as  E^  then  o  ,,  m  (l-q)o.  We  note  in  particular  that 

<  0  if  q>l,  which  can  occur  for  certain  ranges  of  E/p  in  CH^  or  in  Ar 
doped  with  CjFq .  (See  Ref.  3  for  examples  of  collision  frequencies  of  various 
gases.)  In  this  case  the  hose  stability  properties  are  unclear  at  present, 
but  preliminary  investigations  indicate  that  low-frequency  modes  are  damped 
but  high-frequency  modes  may  be  unstable  (the  opposite  of  the  usual 
situation) . 

It  is  also  possible  for  q<0,  e.g.,  the  Spitzer  regime, 
CO.,.  In  this  case  stability  should  be  improved. 

There  are,  of  course,  monopole  effects  due  to  a 
temperature  dependent  chemistry  model  which  affect  the  hose. 
The  higher  temperatures  in  the  nose  region  increase  the 
collision  frequency.  This  will  depress  the  conductivity  and 
cause  more  hose  growth.  On  the  other  hand  the  opposite 
effect  happens  in  the  tail  where  it  is  cooler. 

The  radial  profile  of  the  temperature  is  monotonically 
decreasing  in  most  cases.  In  air,  this  effectively  broadens 
the  conductivity  which  is  stabilizing  for  the  hose.  This 
effect  is  pronounced  in  density  channels  where  the  lower 
density  on  axis  results  in  high  temperatures. 

2 

4.  Spread  Mass  Relation.  Using  the  spread  mass  model  of 
Lee,  one  can  derive  a  dispersion  relation  equation  for  the 
hose  when  the  collision  frequency  varies  as  a  power  of  E. 

Let  v*E^.  If  we  assume  (t)  direct  beam  ionization  is 
the  only  source  of  electrons,  (2)  the  return  current 
fraction  f-I  /I&  where  I  is  the  plasma  current  is  fixed  and 
(3)  the  radial  profile  of  q  and  are  Bennett  with  the 

same  radius,  then  the  electric  field  is  determined, 

■  •to  A 


where  k  is  a  constant. 


A--(dA/dr)D  and  J-- 


The  spread  mass  model  assumes  A--(dA/dr)D  and  J-- 
(dJ/dr )D( 1 +G) -where  all  the  ;  dependence  is  in  0  and  G  is  a 
function  of  u  and  Q* .  where  Q  is  the  z-frequency  of  the 
disturbance  and  Q.  B0ls  the  on-axis  betatron  frequency. 
Making  these  substitutions  results  in 
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where 


A 


q 


th«  ease  q-0  has  b<^n  solved  in  Ref. 


4. 


5.  Coda  Results.  *  In  this  aaetion  we  praaant  VIPER  code 
results  illustrating  the  theory.  The  collision  frequency 
used  in  VIPER  is  the  sua  of  the  electron-neutral  collision 
frequency  of  Fig.  1  and  the  Spitzer  term.  The  electron 
teaperature  in  eV  is  given  by  an  E/P  nodal; 

T  -  O.KE/P)0,8  +  T 

where  E/p  is  in  volts/cn/torr  and  T°  is  a  constant  nodeling 
the  direct  beaa  heating.  The  other  Sera  was  derived  by  All" 
froa  ionization  rates. 

Figure  2  shows  the  aaziaua  hose  diaplaceaent  (saturation 
had  not  occurred.  by  the  tine  of  the  graph)  as  a  function  of 
C  for  a  10  kA,  5  aa  bean.  The  three  curves  correspond  to 
three  VIPER  runs.  In  A  the  full  dipole  conductivity  terns 
are  retained.  In  B  the  dipole  conductivity  is  set  equal  to 
aero.  In  C  only  the  contribution  to  £  due  to  S  is  retained. 
As  seen  in  the  figure  this  part  of  the  dipole  conductivity 
is  stabilizing.  This  is  explained  in  Ref.  4.  The  inclusion 
of  the  w  tera,  Case  A,  is  seen  to  be  destabilizing  in 
accordance  with  theory. 

Finally,  Fig.  3  shows  the  on-axis  conductivity,  <j, 
electric  field  and  fraction  s  ff/e  as  a  function  of  c.  For 
this  beaa  the  conductivity  ®1S  reduced  by  nearly  50%  in 
regions  where  the  hose  is  dangerous. 
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DIAGNOSTIC  DEVICES 


M.  Raleigh,  R.E.  Pectaacak,  S.  Eauver,  D.P.  Murphy,  and  J.R.  Greig 

Saral  Research  Laboratory 
Experimental  Plasma  Physics  Branch 
Plasm  Physics  Division 
Washington,  DC  20375 


INTRODUCTION 

Daring  the  course  of  our  CPB  propagation  experiments  ve  have  found  it 
necessary  to  design  a  number  of  rather  simple,  hut  none  the  less  effective  and 
useful,  diagnostic  systems*  In  this  manuscript  ve  describe  first  a  single 
frame,  gated,  image  intensified  camera  and  then  some  very  fast  B-dot  loops 
that  are  hack  terminated  so  that  several  of  them  can  ganged  together  to 
produce  super-fast  Ibgovski  belts  or  devices  sensitive  to  the  position  of  the 
current  centroid.  These  are  not  new  inventions  but  rather  conveniently 
packaged  versions  of  sell  known  devices  for  which  working  drawings  are 
available. 


Thu  GATED  MICROCHANNEL  PLATE  REFLEX  CAMERA 

The  original  purpose  of  this  camera  was  to  make  it  possible  to  photograph 
a  typical  REB/ atmosphere  interaction  in  the  presence  of  a  preformed  luminous 
channel.  The  camera  (Figure  1)  consists  of  a  standard  Haaselblad  500/C  camera 
with  a  Iblaroid  film  back  which  has  been  modified  so  that  a  proximity  focussed 
MicroChannel  plate  intensifler  tube  (ITT  PUlll)  is  positioned  between  the  film 
plane  of k the  camera  and  the  Iblaroid  film.  This  intensified  tube  has  been 
gated  to  exposure  times  of  a  few  nanoseconds1  and  as  mounted  in  this  camera 
should  be  capable  of  exposure  times  as  low  as  -  10  ns.  So  far  we  hare 
demonstrated  exposure  times  down  to  -  25  ns.  The  frame  size  is  a  circle 
-  1  cm  diameter  and  the  optical  gain  of  the  system  is  ~  1000. 

The  image  intensifler  package  consists  of  three  parts:  a  metal  adapter 
plate  that  mates  with  a  Hasselblad  500/C,  an  aluminum  housing  that  contains 
the  intensifler  tube,  and  a  standard  plastic  Iblaroid  film  pack  holder.  The 
adapter  plate  and  the  film  pack  holder  are  parts  from  a  standard  500/C 
Iblaroid  back.  In  its  operating  position,  the  image  tube  photocathode  is 
located  at  the  image  plane  of  the  lens  and  the  output  fiber-optic  plate  is 
pushing  gently  against  the  Iblaroid  film.  The  film  pull-tabs  are  covered  by  a 
lever.  Pulling  the  lever  to  uncover  the  pull-tabs  also  moves  the  intensifler 
tube  away  from  the  film  allowing  the  film  to  be  pulled  without  wiping  on  the 
fiber-optic  plate.  A  scene  to  be  photographed  through  the  image  intensifler 
is  focussed  through  the  reflex  optics  of  the  camera.  No  continuous  operation 
'focus*  mode  is  necessary  for  the  image  intensifler.  Further,  by  replacing 
the  image  intensifler  package  with  a  Iblaroid  film  back,  the  scene  can  be 
photographed  directly  for  alignment  or  reference  purposes.  Lastly,  this  image 
Intensifler  system  has  the  advantage  of  a  complete  camera  system:  a  variety 
of  commercially  available  accessories  and  a  large  variety  of  lenses  each  with 
its  own  mechanical  shutter  and  iris. 


8D4PLE  B-DOT  LOOPS  AHD  R0G0W5KI  BELTS 


Tso  specific  B-dot  loops  have  Been  used  at  URL.  (Figures  2b  and  3b).  The 
larger  loop  (Figure  2b)  has  an  area  of  -  5  car  and  Is  a  single  turn  loop 
alaoet  totally  enclosed  In  3/l£  Inch  soft  copper  pipe.  Ihe  loop  is  "back 
terminated"  as  shown  with  a  50  Q  resistor  to  prevent  reflections  vhen  several 
(typically  eight)  of  theee  loops,  situated  around  the  circumference  of  a 
circle,  are  used  with  the  sunning  network  shown  in  Figure  2a  to  fora  a 
Bogovski  belt.  The  rlsetiae  of  these  loops  and  of  the  60  ca  4 leant er  Rogovski 
belt  containing  eight  loops  is  ~  T50  ps. 

The  snaller  loop  (Figure  3b)  has  an  area  of  only  -  0.5  ca2,  is  a  single 
turn  loop,  and  is  embedded  as  indicated  in  a  block  of  »"»•*"»-_  Again  the 
loop  is  "back  terminated"  with  a  50  0  resistor  so  that  the  outputs  for  several 
loops  can  be  summed  to  fora  a  Rogovski  belt.  In  this  case  four  or  eight  such 
loops  are  embedded  in  a  circular  groove  cut  in  an  aluminum,  end  flange  (Figure 
3a)  through  which  the  KEB  is  made  to  propagate.  These  loops  and  their 
Ibgowski  belts  have  a  measured  risetime  of  ~  250  ps.  Four  loop  belts  on  a  5 
ca  disaster  and  eight  loop  belts  on  a  25  ca  diameter  have  been  used. 

Vhen  theee  Rogovski  belts  are  used  in  conjunction  with  our  fast 
integrators  (ERL  Meaorandua  Report  Bo.  4939  (1902))  which  have  a  rlsetiae 
of  ~  300  pe,  net  current  measurements  can  be  aade  with  time  resolution 
compatible  with  that  achieved  with  a  Tektronix  710 U  oscilloscope,  l.e. ,  a 
total  system  rlsetiae  of  ~  500  ps. 
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THE  DIAGNOSTICS  0."  HIGH  CURRENT  GENERATED  AIR  PLASMAS 
A.  W.  All,  Naval  Research  Laboratory,  Washington,  D.  C. 


1.  INTRODUCTION 

The  propagation  of  a  high  current  electron  beam  In  air  Is  accompanied  by 
the  air  plasma  and  the  heating  of  air.  Thus  a  channel  Is  generated  with 
numerous  properties  that  control  the  propagation  and  the  guidance  of  the  beam. 
These  properties  are  time  and  space  dependent  and  depend:  on  various 
properties  of  the  beam  current.  In  general,  however,  one  can  divide  the 
channel  characteristics  Into  four  distinct  time  scales;  which  are,  the  early 
time,  the  saturated  region,  the  beam  turn  off  region  and  the  late  time.  The 
properties  of  the  plasma  and  the  channel,  such  as  the  conductivity,  etc.,  are 
important  parameters  which  must  be  determined  over  the  entire  time  scale  for  a 
comprehensive  understanding  of  the  channel  physics. 

In  this  paper  we  discuss  the  measurements  of  several  spectroscopic 
quantities  which  will  provide  quantitative  answers  to  the  channel  properties 
of  interest  for  a  high  temperature  plasma. 

2.  CHANNEL  PLASMA  PROPERTIES 

The  following  quantities  are  considered  as  Important  channel  properties 
to  be  determined:  (a)  The  conductivity,  (b)  The  electron  density,  (c)  The 
electron  temperature,  (d)  The  degrees  of  Ionization  and  dissociation,  (e)  The 
heavy  particle  temperature  and  (f)  The  vibrational  temperature. 

These  quantities  can  be  obtained  from  certain  spectroscopic  measurements 
discussed  in  the  next  section. 

3.  SPECTROSCOPIC  MEASUREMENTS 

We  propose  the  time  resolved  and  space  resolved  meansurements  of  the 
following  quantities: 
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(a)  The  intensity  at  3371A  and  3914A  from  N2  and  N2+ 

(b)  The  Intensity  and  the  Stark  width  of  5330A  and  6157A  from  oxygen  atom 

(c)  The  Intensity  and  the  Stark  widths  of  4222A  and  4256A  from  nitrogen 
atom 

(d)  The  intensity  and  the  Stark  width  of  one  or  more  of  these  lines  at 

4913A,  4341A,  4111A,  4074A  and  3966A  from  oxygen  ion  (0+) 

(e)  The  Intensity  and  the  Stark  widths  of  one  or  more  of  these  lines  at 

4613A,  4643A,  5006A,  4447A  and  3995A  from  nitrogen  ion  (N+) 

(f)  And  the  intensity  of  a  20A  band  continuum  emission  at  4834A 

These  Intensities  and  widths  provide  abundant  information  on  the  electron 
density,  electron  temperature,  atomic  and  molecular  densities  and  the 
densities  of  their  respective  ions. 

4.  CALCULATIONS  OF  THE  MEASURED  QUANTITIES 

The  Stark  widths  for  the  lines  suggested  above  have  been  calculated  and 
measured  for  specific  electron  densities  and  temperatures.  Their  data  can  be 
Interpolated  or  extrapolated  for  any  desired  electron  density. 

The  lines,  bands  and  the  continuum  Intensities  can  be  calculated  using 
detailed  ablnltlo  electron  energy  deposition,  air  chemistry  and  emission 
codes.  Some  typical  continuum  calculations  are  presented  for  several  plasma 
tenperatures  In  Figures  1,  2  and  3  where  the  free-free  (electron -Ion),  the 
free-bound,  and  the  free-free  (electron -neutral)  emissions  are  shown.  In 
Figure  3  Indications  are  given  for  band  emissions  from  NO,  H20  and  C02  in  the 
Infrared  region  which  will  be  superimposed  over  the  calculated  continuum. 
These  band  emissions  can  be  calculated  by  the  knowledge  of  the  species 
densities,  the  vibrational  and  the  gas  temperatures  In  the  channel  provided  by 
air  chemistry  models. 


S.  CONCLUSIONS 


With  a  redundant  Information,  from  the  spectroscopic  measurements,  on  the 
electron  density,  electron  temperature  and  species  densities,  a  chanllenging 
opportunity  arises  for  the  code-experimental  data  analysis  and  interaction. 
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ABSTRACT 

We  describe  an  experiment  in  which  an  Intense  REB  has  been  allowed  to  interact 
with  reduced  density  channels  in  ■iwnnls  at  an  ambient  pressure  of  40  Tbrr.  Ammonia 
was  chosen  for  these  experiments  because  of  its  ability  to  absorb  the  energy  from  a 
pulsed  COg  laser.  This  permitted  the  formation  of  reduced  density  channels  in  an 
otherwise  uniform  atmosphere  of  mil  a.  At  this  pressure  and  with  these  particular 
RIB  characteristics,  avalanche  ionisation  was  not  a  major  source  of  plasma  conductivity 
in  the  ambient  gas.  But  in  the  reduced-density  channel  and  using  high  inductance 
boundary  conditions  avalanche  ionisation  was  the  major  source  of  plasma  conductivity, 
easily  exceeding  beam  induced  conductivity  (i.e. ,  by  direct  particle  collisions).  Then 
the  presence  of  the  reduced-density  channel  significantly  Influenced  REB  propagation, 
Plaaaa  return  currents  were  markedly  increased  and  confined  to  the  channel  with  the 
result  that  the  RIB  was  always  expelled/repelled  from  the  channel, 

APPARATUS 

These  propagation  experiments  were  performed  in  a  large  test  chamber  made  of 
fiberglass.  The  chamber  was  cylindrical:  1,5  a  long  and  0-6  m  in  disinter  (Figure 
1),  The  chafer  could  be  fitted  with  a  brass  conducting  liner  to  ehame  from  the 
initial  nonconducting  mall  boundary  conditions  to  conducting  wall  boundary 
conditions.  The  COg  laser  beam  entered  the  test  ehanber  through  a  salt  (BaCl)  window 
on  the  end  of  a  2.5  a  long  extension  tube.  Further  details  are  available  in  IRL 
Memorandum  Report  #$261  (1984). 

CHAHEL  AID  REB  CHARACTERISTICS 

%  create  reduced  density  channels  in  the  ammonia  gas  a  10  a  focal  length,  salt 
lens  focused  the  COg  laser  been  through  a  circular  aperture  on  the  salt  window  that 
clipped  the  laser  beam  to  a  2  cm  radius  spot  where  the  REB  entered  the  test  chamber. 

The  heated  gas  expanded  to  fora  a  simple  channel  of  radius  ~  2.5  cm  (Figure  2). 

Because  the  effects  of  absorption  and  focusing  were  nearly  matched,  the  laser 
energy  deposited  in  the  smmnnla  gas  was  fairly  uniform  at  ~  25  mJ/caJ along  the  length 
of  the  channel.  This  amount  of  energy  (~  0.1  eF/molecule)  while  several  times  the 
ambient  energy  (~  0.025  eV/nolecule)  was  negligible  compared  to  the  dissociation 
potential  (4.3  eV)  and  the  ionisation  potential  (10.15  eV)  of  SMonla. 

The  vibrational-translational  relaxation  of  the  excited  ammonia  molecule  is  very 
f**t«  ~  100  ns  at  a  pressure  of  40  Tbrr.  Thus  we  expected  channel  formation  to  result 
from  effectively  Instantaneous  tberaalisation  followed  by  relaxation,  on  hydrodynamic 
time  scales,  to  pressure  equilibrium.  We  did  not  expect  the  laser  to  generate  any 
conductivity  in  the  channels.3  These  aspects  of.  the  channel  history  were  confirmed 
using  double  exposure  holographic  interferometry4  and  an  RF  bridge*  respectively. 
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The-Pulserad  310  generator  produced  an  REB  of  approximately  19  kA  peak  current  at 
1.1  MV  peak  voltage  with  a  pulse  *tdth  of  26  ns  FVHM.  The  REB  was  transported  to  the 
test  chafer  through  a  drift  tube  jntainlng  a  conical  return  current  conductor  and 
filled  with  dry  nitrogen  at  20  Tbr.  It  entered  the  test  ehaaiber  through  a  2.3  cm 
aperture.  The  peak  REB  current  entering  the  test  chamber  was  ~  8  kA  and  it  had  the 
current  density  profile  shown  in  Figure  2. 

EXPERIMENTAL  RESULTS 


The  REB 
contitlons: 


was  Injected  into  the  test  chamber  under  four  distinct  experimental 


(1)  conducting  liner  installed  -  no  channel  present 
(ll)  conducting  liner  installed  -  channel  present 
(lii)  conducting  liner  removed  -  no  channel  present 
(iv)  conducting  liner  removed  -  channel  present. 


The  time  histories  for  the  measured  REB  and  net  currents  along  with  the  calculated 
plasma  current  are  shown  in  Figure  3  for  case  (l).  Bote  that  the  plasma  current  was 
negative  during  mich  of  the  REB  pulse.  In  all  four  cases  the  REB  failed  to  propagate 
the  length  of  the  test  chamber  and  deviation  from  its  initial  trajectory  increased  with 
each  case:  (1)  through  (iv).  Ve  were  able  to  correlate  the  magnitude  of  the  deviation 
with  a  parameter  ve  call  f*  where 

<V_. 


Case 

REB  Behavior 

1 

**  *3 

Moderate  hosing 

ii 

-  .3 

Moderate  hosing 

ill 

-  .5 

Violent  hosing 

iv 

-  .7 

REB  dispersed 

AIALX8I8 

The  conductivity  of  the  gas  in  the  chamber  is  a  function  of  the  electron  density 
and  the  gee  density,  e  ■  e(n  /n).  The  growth  of  conductivity  was  governed  by  the 
production  rate  of  electrone*generated  by  the  passage  of  the  REB  through  the  gas.  Four 
contributions  to  the  rate  equation  were  direct  eolllslonal  ionisation,  avalanche 
ionisation,  dissociative  attachment  and  dissociative  recombination. 
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four  terms  are  illicit  and/or  explicit  functions  of  the  gas  density,  n,  the 
electron  density,  n#,  the  electric  field,  E,  and  the  REB  current  density,  Jg.  The 
avalanche  ionisation  coefficient,  e,  is  in  particular  a  very  strong  function  of  the 
density  normalised  electric  field,  E/n.  Figure  b  is  a  plot  of  the  density  normalised 
avalanche  coefficient,  s/n,  as  a  function  of  E/n  for  both  air  and  ammonia.3*6  The 
avalanche  coefficient  increases  if  either  the  gas  density  drops,  as  in  a  channel,  or 
the  electric  field  increases.  Our  REB,  when  injected  into  the  lined  test  chamber  with 
bO  Tbrr  eautnia  present  produced  an  initial  density  nornalised  electric  field 
of  ~  200  Td.  This  seas  REB  injected  into  a  reduced  density  channel  would  develop  a 
larger  density  normalised  electric  field  of  -  800  Td. 
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Baiii(  in  tha  limit  where  the  risetime  of  the  REB  was  longer  than  any  transit  time 
for  electromagnet ic  waves  within  the  test  system,  the  electric  field  generated  by  the 
HEB  was  E  *  L  I|  where  I.  is  the  net  current  and  L  is  the  inductance  per  unit  length 
scan  by  tha  REB.  Whan  tha  tast  chamber  was  lined  with  a  conducting  screen  this  luiqped 
circuit  inductance  was  estimated  to  be  that  for  two  concentric  coaxial  cylinders  of 
radii  a  ■  1.2  cm  and  b  *  30  cm. 


L  "  "  a  la(f) 


•5l<0  (H/m) 


Dm  tast  chamber  sat  in  a  metal  walled  room  of  3  x  3  m  cross  section.  When  the  screen 
liner  was  removed  from  tha  teat  chamber  the  outer  conducting  boundary  radius,  b, 
defaulted  to  that  of  the  metal  experimental  room,  b  ~  ISO  cm.  The  estimated  ltuped 
circuit  inductance  for  the  nonconducting  test  chamber  was  L  , 

QC 


Since 
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the  density  normalised  electric  field,  E/n, 
tug  liner  was  removed  from  the  test  chamber. 


was  larger  when  the 


We  developed  a  counter  model  of  the  REB/plasma  interaction  from  a  simple  linked 
parameter  circuit  model.  The  model  used  the  measured  time  history  of  the  REB,  the 
channel  parameters  and  the  boundry  conditions  to  calculate  the  plasma  and  net 
currents.  The  model  predicted  values  of  f*  very  close  to  those  actually  measured  in 
the  experiment  for  sets  of  parameters  that  matched  the  true  experimental  conditions. 

The  effect  of  changing  tha  experimental  conditions  warn  to  alter  the  relative 
importance  of  the  terms  in  the  rate  agnation.  For  enable,  in  case  (iv),  the  avalanche 
ionisation  rate  in  tha  channel  was  increased  by  approximately  two  orders  of  magnitude 
over  ease  (l)  sad  dnmi noted  the  other  terms  in  the  electron  rate  equation.  The  model 
predieted  a  large  increase  la  the  channel  conductivity  and  a  very  localised  plasma 
return  currant,  exactly  as  observed  in  the  experiment. 

COKUBIO* 

We  have  injected  an  Intense  REB  into  a  test  chamber  filled  with  bO  Tbrr  ammonia 
gas  and  observed  the  effect  of  boundary  conditions  and  reduced  density  channels  on  the 
propagation  of  the  REB.  Ws  have  also  shown  that  a  simple  co^xiter  model  of  the 
interaction  accurately  forecasts  the  most  important  macroscopic  parameters  of  the 
interaction. 

tftider  the  conditions  of  an  enhanced  density  normalised  E  field,  such  as  in  case 
(iv)  where  f*  ~  .7,  propagation  of  the  REB  was  co^letely  disrupted  with  the  REB 
striking  the  chamber  walls  within  30  cm  of  injection.  When  the  experimental  conditions 
were  adjusted  so  the  avalanche  term  was  less  dominant  (f*  ■  0.3  to  .0.5)  the  REB  was 
again  ejected  from  the  channel  and  was  hose  unstable  but  did  not  break  apart.  At  no 
time  did  the  simple  density  channels  attract  or  guide  the  REB;  their  effect  was  always 
disruptive. 


This  work  was  supported  by  the  Defense  Advanced  Research  frojects  Agency  and  the 
Office  of  Rkval  Research. 
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